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Abstract 

The fluctuations of tlie largest fragment cliarge of a partition and of the 
charge asymmetries of the two or three largest fragments in spectator decays 
O ! following 19'^Au + i^'^Au colhsions at 1000 MeV per nucleon are investigated. 

■ The observed bimodal distributions at specific values of the sorting variable 

Zbound exhibit features known from percolation theory where they appear as 
finite-size eff^ects. The underlying configurational fluctuations seem generic 
for fragmentation processes in small systems. 



X 



^3 



o 



1 Introduction 



Double-humped event distributions have received particular interest recently be- 
CJ ' cause their observation might indicate bimodality which is one of the signals ex- 

pected for a first-order phase transition in finite systems [1-3]. Bimodality occurs 
when non-negligible surface interactions at the phase boundary lead to a convex 
entropy function in the transition region. Canonical sampling near the transition 
I temperature will then produce two distinct event classes which differ with respect 

to the order parameter of the transition. 

^ ■ Bimodality has, e.g., been observed in solid-liquid transitions of clusters of Na 

O 



atoms In these experiments, the clusters were thermalized in a heat bath of 
helium gas and excited with photons from a laser beam. Corresponding experi- 
I ments in nuclear fragmentation face the difficulty that the temperature cannot be 

^ I predetermined and that a canonical sampling can thus not be performed. There is, 

furthermore, the possibility of impact-parameter mixing, meaning that fluctuations 
^ . of the variables used for sorting will unavoidably lead to flnite distributions with 

respect to other observables or event characteristics (e.g., impact parameter) even 
in narrowly selected event samples. It is, nevertheless, of interest to study the origin 
and the meaning of bimodal event distributions and the conditions under which they 
occur. For example, Lopez et al. using the HIPSE event generator have recently 
pointed to the role of angular momentum in producing bimodal distributions as a 
result of instabilities of nuclei with high spin [5]. Pichon et al. [3] have shown that 
the two bumps of the bimodal distributions observed in ^^^Au + ^^^Au fragmenta- 
tions at 60 to 100 MeV per nucleon correspond to different scaling properties of the 
distributions of the largest fragment charge (A-scaling [0]). 

In this work, the fluctuations of the largest fragment charge of a partition and 
of the charge asymmetries of the two or three largest fragments from the decay of 
excited projectile spectators in ^^''Au + ^^''Au collisions at 1000 MeV per nucleon 
are investigated. The data have been collected in experiments performed with the 
ALADIN spectrometer at GSI. The observed bimodal distributions at speciflc values 
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Figure 1: Distribution of Z^^.^ ver- 
sus Abound for projectile fragments from 
Au on Au at 1000 MeV per nu- 
cleon [15]. Conventional fission events 
are removed. The shadings follow a log- 
arithmic scale. 
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Figure 2: Bond percolation: mean rel- 
ative magnitude of the largest cluster 
as a function of the bond probability 
for cubic lattices of Cq = 81 and 27000 
sites. The critical bond probability in 
the infinite sytem is = 0.249. 



of the sorting variable Abound exhibit features known from percolation theory where 
they appear as finite-size effects. Percolation on a large lattice exhibits signatures of 
a second-order phase transition |7j . The observed similarities thus raise the question 
whether bimodality in fragmentation reactions may be used to infer the order of the 
phase transition in the nuclear case. 

Bimodality and its significance for the interpretation of fragmentation data and 
their relation to the nuclear liquid-gas phase transition have been frequently dis- 
cussed at previous conferences of this series [8-12]. At this year's conference, new 
results obtained by the INDRA collaboration are presented by E. Bonnet [T3] . 

2 Experimental results 

The data used for the present analysis were obtained by the ALADIN collaboration 
in measurements using ^^''Au projectiles of 1000 MeV per nucleon delivered by the 
heavy-ion synchrotron SIS at GSI [T11|T3]. The ALADIN spectrometer was used to 
detect and identify the products of the projectile-spectator decay following collisions 
with ^^^Au target nuclei. 

The sorting variable .Abound is defined as the sum of the atomic numbers Zi 
of all projectile fragments with Zi > 2. It reflects the variation of the charge of 
the primary spectator system and is monotonically correlated with the impact pa- 
rameter of the reaction [16]. The evolution of the dominant reaction processes is 
illustrated in Fig. [1] which shows the correlation of the largest atomic number Zmax 
observed in a partition with Abound- Large values of Abound correspond to low ex- 
citation energies, at which the decay changes its character from evaporation-like 
processes ~ Abound) to multifragmentatiou ("rise" of multifragmentation) 



while small values correspond to reaction channels with high excitation energies and 
disintegrations into predominantly very light clusters ("fall" of multifragmentation, 

•^max -Abound ) • 

Besides the evolution of the mean and of the fluctuations of ^max (alternatively 
denoted by in the following) also those of two-fragment and three-fragment asym- 
metries are of interest and characterize the dominant transition of the reaction mech- 
anism pTT]. It is found that the ratios {Z2/Z1) and (Z3/Z2) both approach 0.6 
at small .Abound {Z2 and Z^ are the second and third largest atomic number Z of a 
partition). Consequently, the charge difference Z\ — Z2 — Z^, or the corresponding 
asymmetry after normalizing with respect to the system charge Zq, will approach 
zero at small 2'bound while it is close to Zq (the asymmetry close to 1) at large 2'bound- 
In the transition region, the fluctuations of these observables are large [17J and the 
distributions are bimodal, i.e. they exhibit a two- hump structure (Fig.[3l). Note that 
(-^max) drops most rapidly in the bin 53 < Abound < 57 at which this bimodality is 
most strongly pronounced (Fig. [T]). 
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Figure 3: Distributions of the atomic number Zi of the largest fragment of a partition 
and of the differences Zi — Z2 and Z\ — Z^ — ^3 for events with 53 < .Abound < 57 
from the fragmentation of ^^^Au projectiles at 1000 MeV per nucleon. Note the 
offsets by 200 and 400 counts of the difference distributions. 



3 Largest fragment as order parameter 

In the search for an experimentally accessible order parameter of the nuclear liquid- 
gas phase transition, as observed in multifragmentation reactions, the magnitude of 
the largest fragment of the partition has appeared as a promising choice. It may 
be identified with the part of the system in the liquid phase, and it is correlated 
with the mean density which is the natural order parameter of a liquid-gas phase 
transition. Observables correlated with it, as e.g. the differences and asymmetries 
discussed above, may similarly serve as order parameters. 



Statistical model calculations for nuclear multifragmentation show that the dis- 
appearance of the dominating fragment is associated with a maximum of the heat 
capacity which is the more strongly pronounced the larger the system [18]. For 
Aq = 150, the system mass for Zbound ~ 50 [H], the predicted specific heat distribu- 
tion is rather wide with a maximum at T 6.3 MeV. This transition temperature, 
or boiling temperature according to the authors of [TH], is comparable with values of 
the double-isotope temperature TneLi as measured for the present reaction [HI [19] 
and for similar systems [20]. The good description of the charge correlations and 
charge asymmetries characterizing the partitioning of the system, including their 
variances, with statistical multifragmentation models provides further evidence for 
the first-order nature of the transition [THl 121] • Bimodality is predicted for canonical 
ensembles [22] [23] . 

The disappearance of the largest cluster, with the variation of a suitable control 
parameter, has been identified as a prominent signal also in fragmentations of other 
systems as, e.g., atomic hydrogen clusters ^], and the extension of the largest 
cluster is an order parameter in percolation theory [7] . On finite percolation lattices, 
the disappearance of a dominant largest cluster proceeds rather smoothly and with 
obvious similarity to the nuclear experiment (Figs. [T][2|) . 
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Figure 4: Examples of the distributions of the 3-fragment difference Ci — C2 — C3, 
normalized with respect to the lattice size Cq = L^, as obtained with cubic bond 
percolation for the cases Cq = 125 sites and pb = 0.328 (top) and Cq = 27000 sites 
and Pb = 0.259 (bottom). 



4 Cubic bond percolation 



Percolation models have been quite successfully used for describing the properties 
of nuclear fragmentation [17,24-28] including the apparent critical behaviour. For 
the present purpose, calculations with a cubic-bond-percolation model have been 
performed with various lattice sizes. The critical bond parameter for this type 
of lattice is = 0.249 [71 ^U\. For large lattices, in the limit of infinity, a sharp 
transition with the sudden appearance of an extended percolating cluster is observed 
for this value of the probability that a bond exists between neighbouring sites. For 
finite lattices, the transition is smooth (Fig. [2]) and, for a lattice of 81 sites (obtained 
by smoothing the corners and edges of a 5^ lattice for simulating the 79 charges of a 
Au nucleus) it is very similar to what is observed in the nuclear experiment (Fig. [1]). 
For specific values of the bond parameter in the transition region, the distributions of 
the 3-cluster asymmetry Ci — C2 — C3 (the cluster sizes Cj are ordered in magnitude) 
exhibit two bumps (Fig. H]). For the smaller lattice of 125 sites, the distribution 
extends over a major part of the interval [0,1] that is accessible after normalization 
with respect to the number of sites cq = L^. Also this feature is reminiscent of the 
result obtained for the ^^''Au fragmentation (Fig. [3]). 
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Figure 5: Bond parameter Pbimod for which the bimodal patterns appear most clearly 
(top) and difference Pbimod — p'^ (bottom) as a function of the lattice size L. The 
full hues represent the results of three-parameter (top) and two-parameter (bottom, 
with fixed) power-law fits according to Eq. [TJ The dashed line indicates the 
location of the critical bond parameter p'^ = 0.249 for the infinite system [TJ [SO] . 




For the larger lattice, the distribution is still double-humped but becomes much 
narrower. The bond parameter Pbimod at which the bimodal structure is most pro- 
nounced is smaller and much closer to the critical value. Calculations performed for 
various lattice sizes and samples of up to 200000 events show that this variation is 
systematic and confirm that the law of finite-size scaling is obeyed by Pbimod- A 
power law fit according to the expression 

PUmod-pT = C-L-^^'' (1) 

shows that the critical bond parameter for the infinite lattice is indeed approached 
by Pbimod (Fig-E]). A two-parameter fit with a fixed p'^ = 0.249 yields u = 0.88±0.01 
in agreement with the known value u = 0.88 of the critical exponent describing the 
divergence of the correlation length 

cubic bond percolation 
I — I — ' — I — ' — I — I — I — I — I — I — I — I 




0.1 0.2 0.3 0.^ 0.5 0.6 



Pb 



Figure 6: Root-mean-square widths of the Ci (full line, blue), Ci — C2 (dashed), and 
ci — C2 — C3 (full line, red) distributions as a function of the bond probability pb- 

Finite-size scaling identifies the observed phenomenon as originating from order- 
parameter fluctuations near the percolation phase transition. The same law with 
the same exponent is also valid for the locations of the maxima of the slopes of 
the (ci) vs. Pb transition which practically coincide with Pbimod (cf- Fig. [5]), or 
for the widths of the transition region [7j. The fluctuations of the largest cluster 
size cause corresponding fluctuations of the differences and asymmetries (Fig. [6]), a 
property that identically appears in the fluctuation widths of the charge differences 
or asymmetries observed in the fragmentation of ^^^Au (Fig. [7]). The existence of 
two bumps in the event distributions appears as a generic feature of fragmentation 
processes, including that modeled with percolation. The mere observation of this 
phenomenon can thus not be considered as providing evidence for a first-order phase 
transition. 
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Figure 7: Root-mean-square widths of the Z\ (full line, blue), Z\ — (dashed), 
and Z\ — Z2 — Z^ (full line, red) distributions as a function of Abound- Normalization 
with respect to 2'bound will shift the maxima to Abound ~ 55, i.e. into the center of 
the transition region (cf. Fig. [T]). 

5 Reaction scenarios 

In classical molecular dynamics, maximum size fluctuations define a critical per- 
colation line (Kertesz line), or a critical percolation band in finite systems, in the 
temperature-density phase diagram [32] • The Kertesz line, known from studies of 
the lattice-gas model [HHllHlj, extends from the thermodynamical critical point into 
the supercritical region of higher density and temperature and is considered generic 
for simple fluids. Its identification requires appropriate algorithms for the recogni- 
tion of clusters in the dense medium. Equilibrium cluster-size distributions along 
the critical band exhibit a power-law behaviour and bimodality [321 

It is a particular characteristic of the classical-molecular-dynamics model that 
the distributions of so-defined clusters do not significantly change as the systems are 
allowed to expand freely to a breakup point beyond which clusters can be recognized 
in configuration space [3^. Their properties acquired by originating from a phase 
space location in the critical region will be reflected in the asymptotic distribu- 
tions. The reaction scenario suggested by these calculations thus links the observed 
percolation-like phenomena to a truly critical behaviour of large systems. The ap- 
plicability of the model to nuclear fragmentation can be tested by searching for 
predicted non-equilibrium phenomena at breakup [32]. One of them, a considerable 
difference between the internal temperatures of the emerging fragments and that 
of the environment, is also a result of quantum-molecular dynamics (QMD [36j). 
The recent analysis of multifragmentation following ^^''Au -|- ^^''Au collisions in the 
energy range 60A to 150A MeV with this model has, in particular, also shown that 
bimodality is observed and that the experimental asymmetry distributions of the 



largest fragments are reproduced rather well with QMD [36j. 

On the other hand, phenomena resembling critical behaviour as it appears in 
large systems are also observed for equilibrium distributions of small systems gen- 
erated within their coexistence zones. For the lattice-gas model, it has been rather 
generally shown that the observation of scaling inside the coexistence zone is compat- 
ible with a first-order phase transition because of finite size-effects • The scaling 
will disappear in large systems. Similar conclusions were reported by the authors 
of [38]. When the Statistical Multifragmentation Model was used to describe the 
fragmentation of relativistic ^^''Au projectiles, the experimentally observed power- 
law Z and bimodal .^max distributions in the transition region have been reproduced 
with conditions below the critical point of this model [21]. These phenomena thus 
seem to appear naturally when viewing phase transitions in small systems through 
their partitioning into fragments. 

An interesting experimental observation is the coincidence of several signals con- 
sidered indicative of a phase transition in fragmentation data [Ml HO]. Besides bi- 
modality, this includes universal fluctuations (A-scaling [6J) of the size of the largest 
fragment and the kinetic-energy fluctuations which have been associated with neg- 
ative heat capacity [H]. Searching for a common origin, it seems most likely that 
they are all related to the configurational fluctuations [121 US] to be expected in 
fragmentation processes and identified as finite-size effects in percolation. 

6 Summary 

Experimental results regarding the largest fragment charge and the asymmetries of 
the two and three largest fragments in the decays of ^^^Au projectile spectators at 
1000 MeV per nucleon have been presented. The bimodal distributions at specific 
values of the sorting variable 2'bound reflect the size fluctuations of the largest 
fragment in the transition region between the regimes of residue production and 
of multifragmentation. In the reaction scenario suggested by molecular dynamics, 
these configurational fluctuations are related to a critical percolation region in 
the phase diagram which reduces to a critical percolation line (Kertesz line) in 
large systems. In small systems, critical-like phenomena like scaling, power-law 
cluster distributions and bimodality are also exhibited by equilibrium distributions 
generated at locations within the coexistence region. The underlying configurational 
fluctuations, identified as finite-size effects with percolation, thus seem generic for 
fragmentation processes in small systems. 

Stimulating discussions with X. Campi and E. Plagnol are gratefully acknowl- 
edged. 
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